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ABSTRACT 

We present spectra of Eris from the MMT 6.5 meter telescope and Red Chan- 
nel Spectrograph (5700-9800 A; 5 A pix" 1 ) on Mt. Hopkins, AZ, and of Pluto 
from the Steward Observatory 2.3 meter telescope and Boiler and Chivens spec- 
trograph (7100—9400 A; 2 A pix -1 ) on Kitt Peak, AZ. In addition, we present 
laboratory transmission spectra of methane-nitrogen and methane-argon ice mix- 
tures. By anchoring our analysis in metha ne and nitrogen solubilities in on e 



Prokhvatilov &: Yantsevichl (119831 ). 



another as expressed in the phase diagram of 
and comparing methane bands in our Eris and Pluto spectra and methane bands 
in our laboratory spectra of methane and nitrogen ice mixtures, we find Eris' 
bulk methane and nitrogen abundances are ~ 10% and ~ 90% and Pluto's bulk 
methane and nitrogen abundances are ~ 3% and ~ 97%. Such abundances for 
Pluto are consistent with values reported in the literature. It appears that the 
bulk volatile composition of Eris is similar to the bulk volatile composition of 
Pluto. Both objects appear to be dominated by nitrogen ice. Our analysis also 
suggests, unlike previous work reported in the literature, that the methane and 
nitrogen stoichiometry is constant with depth into the surface of Eris. Finally, 
we point out that our Eris spectrum is also consistent with a laboratory ice mix- 
ture consisting of 40% methane and 60% argon. Although we cannot rule out 
an argon rich surface, it seems more likely that nitrogen is the dominant species 
on Eris because the nitrogen ice 2.15 /xm band is seen in spectra of Pluto and 
Triton. 



Subject headings: methods: laboratory — methods: observational — planets and 
satellites: surfaces — techniques: spectroscopic 
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Introduction 



Pluto and Eris are the largest known icy dwarf planets in the Kuiper belt, and the y 
both exhibit methane ice bands in their optical and near-infrared spectra flj 



1993 



Grundy &: Fink 



1996 



Brown. Trujillo. &: Rabinowitz 



2005 



Owen et al. 



Licandro et al 



In addition, Pluto exhibits nitrogen and carbon monoxide ice bands ( 



2006) 



Owen et al. 



1993 



Doute et al. 



19991 ) : however, there are no detections of nitrogen and carbon monoxide bands 



in spectra of Eris. Pluto and Eris have sufficient mass and are sufficiently cold that they 



are able to retain vo 
of the Solar System 



atiles in their atmosphe res, and hence on their surfaces, over the age 



Schaller fc Brown 



20071 ) 



The detection and analysis of methane, carbon monoxide, and nitrogen ice bands 
in Pluto's nea r-infrared sp e ctra has resulted in a number of measurements for Pluto's 



stoichiometry. 



Owen et al. 



( 119931 ) found a methane abun dance of 1.5 %, a ca rbon monoxide 



abundance of 0.5%, and a nitrogen abundance of 98%. 



Doute et al. 



(119991 ) reported a 



methane abundance of 0.5%, a carbon monoxide abundance of 0.1%, and a nitrogen 



abundance over 99%. 



Olkin et al 



(120071 ) found a methane abundance of 0.36%. Nitrogen 



ice dominates methane and carbon monoxide ice on Pluto's surface. 

Because nitrogen and carbon monoxide bands have yet to be detected in spectra 
of Eris, it is more difficult to estimate Eris' stoichiometry. One way to qualitatively 
estimate the methane and nitrogen mixing rati o on E ris is to rely on laboratory studies of 



methane and nitrogen ices. 



Quirico fe Schmitt 



( 119971 ) showed that near-infrared methane 



absorptio n bands shift to sh orter wavelengths when methane is diluted by nitrogen. We 
note that 



Owen et al 



the laboratory work. 



1993 ) saw suc h shift s in their spectra of Pluto four years before 



Brunetto et al. 



( 120081 ) found increasing shifts for methane bands 



with increasing nitrogen content. Eris' methane bands have smaller shifts than Pluto's 
methane bands, suggesting that nitrogen is not as abundant on Eris as it is on Pluto 
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f lBrown. Truiillo. fc Rabinowitzjbooa iDumas et all 120071 Merlin et all 120091 ) 



In addition to using methane band shifts as a proxy for the overall methane- nitrogen 
mixing ratio, it is possible to use the shifts to qualitatively constrain the mixing ratio as a 
function of depth into the surface of Eris. The average penetration depth of a photon at a 
particular wavelength is inversely related to the absorption coefficient at that wavelength, 
e.g. photons corresponding to larger absorption coefficients are absorbed more, preventing 
them from penetrating as deeply into the surface. In other words, stronger bands in the 
spectrum of an icy dwarf planet probe, on average, shallower into the surface than weaker 
bands. So, it is possible to use blue shifts and albedos at maximum absorption of two or 
more methane ice bands to look for a trend in the mixing ratio as a function of depth into 
the surface of Eris. 



Licandro et al.l ( 120061 ) obtained optical spectra of Eris and found evidence of a 
decreasing nitrogen abundance with depth into its surface. Specifically, they found the 
weaker 7296 A band had a shift of 0.1 ± 3 A and the stronger 8897 A band had a shift of 
15 ± 3 A. They speculated that as Eris moved toward aphelion during the last 200 years, 
and as its atmosphere cooled, less volatile methane began to condense out first. As the 
atmosphere cooled further, it became more and more nitrogen rich, until the much more 
volatile nitrogen could condense out on top of the methane rich ice. 



Merlin et al 



( 120091 ) obtained optical and near-infrared spectra of Eris and found 
further evidence for stratification of methane on Eris. Specifically, they found the 7296 A 
and 8897 A bands were shifted by 9 ± 4 A and 12 ± 4 A, respectively. Furthermore, they 
found the stronger near-infrared methane bands, which do not probe a s deep as the optical 



Merlin et al. 



( 120091 ) suggested a pure 



bands, had smaller shifts than the optical bands, 
methane surface layer, on top of a nitrogen diluted methane layer, on top of a pure methane 
layer for Eris. 
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We found similar shifts for five methane bands in the spectrum of Eris ( lAbernathy et al. 



20091 ). For example, we found the 7296 A and 8897 A bands were blue shifted by 7 ± 
1 A and 5 ± 2 A, respectively. Abernathy et al. reported a correlation on the edge of 
statistical significance between blue shift and albedo at maximum absorption for the five 
bands. The correlation suggested an increasing nitr ogen abundance with incre asing depth 



into the surface, i.e. the opposite trend reported by 



Licandro et al. 



fl2006h and 



Merlin et al. 



(120091 ). Abernathy et al. put forth an atmosphere-surface model to explain a way to get an 



increasing nitrogen abundance with depth. 

All of these stratigraphic investigations make the assumption that different methane 
bands should have the same shift for the same stoichiometry. However, laboratory 



experiments show that t he smaller (larger) 



smaller (larger) the shift (jOuirico fc Schmitt 



;he w avelength (wavenumber) of a band the 



19971 ). It is essential to take the different shifts 



for different bands into account before discussing changes in composition with depth. 

Here, we report new optical spectra of Pluto and Eris as well as new transmission 
spectra of laboratory ice mixtures. We use the new data to place quantitative constraints 
on Eris' stoichometry. In addition, we address the issue of whether the composition of Eris 
changes with depth. In sections 2 and 3 below, we describe our astronomical observations 
and laboratory experiments. In section 4, we describe the models we used to derive 
abundances on Pluto and Eris. In section 5, we discuss our results. 



2. Observations 
2.1. Pluto 

We obtained optical spectra of Pluto on the night of 2007 June 21 UT with the Steward 
Observatory 2.3 meter telescope and Boiler and Chivens Spectrograph. We used a 2.5 x 
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240 arc sec entrance slit and a 600 g mm" 1 grating that provided wavelength coverage 
of 7100—9400 A and a dispersion of 2.0 A pixel -1 in first order. There were high, thin 
cirrus clouds and the seeing was ~ 1.5 arc sec throughout the night. Pluto was placed at 
the center of the slit and the telescope was tracked at Pluto's rate. At the start of our 
observations Pluto had an airmass of about 1.6, it transited the meridian at an airmass of 
1.5, and at the end of our observations it had an airmass of about 1.6. 



Since the slit width was 2.5 arc sec and the seeing was ~ 1.5 arc sec, we observed Pluto 
and Charon blended together. Pluto contributed ~ 90% to the total signal. 



Abernathy et al. 



(120091 ). Briefly, we 



We used the calibration techniques described by 
bracketed object spectra with HeNeAr spectra in order to obtain an accurate wavelength 
calibration. Telluric and Fraunhofer lines were removed from the Pluto spectra by dividing 
each 600-second Pluto spectrum by the spectrum of a solar analog star, BS 5996. The 
airmass difference between each Pluto spectrum and its corresponding solar analog spectrum 
was < 0.1. 

We obtained 12 600-second exposures of Pluto. Since we saw no significant difference 
between the individual exposures during the night, we summed them to yield a single 
spectrum of Pluto with an exposure time of two hours. In Figure 1, we present the resulting 
reflectance spectrum of Pluto. The spectrum was normalized to 1.0 between 8200 A and 
8400 Aand moved upward by 0.2 on Figure 1 to facilitate comparison with the spectrum of 
Eris (see below). 

We checked our wavelength calibration by measuring the centroids of ten well-resolved 
sky emission lines across our spectrum and compared them to the VLT high— spectral 



resolution sky line atlas of iHanuschikl (j2003[ ). The average difference and standard deviation 
of the differences for the ten lines were -0.003 A and 0.15 A. So, it appears our sky lines, 
and hence our wavelength measurements, are accurate to ~ 0.2 A. 
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2.2. Eris 

We obtained optical spectra of Eris on the night of 2008 October 3 UT with the MMT 
6.5 meter telescope, Red Channel Spectrograph, and a red sensitive, deep depletion CCD. 
We used a 1 x 180 arc sec entrance slit and a 270 g mm -1 grating that provided wavelength 
coverage of 5700—9800 A and a dispersion of 5.0 A pixel -1 in first order. There were high, 
thin cirrus clouds and the seeing was variable, 0.8 — 1.6 arc sec . Eris was placed at the 
center of the slit and the telescope was tracked at Eris' rate. At the start of our observations 
Eris had an airmass of 1.75, it transited the meridian at an airmass of 1.25, and at the end 
of our observations it had an airmass of 1.89. 



Abernathy et al. 



fl2009[ ). We 



Again, we used the calibration techniques described by 
bracketed object spectra with HeNeAr spectra in order to obtain an accurate wavelength 
calibration. Telluric and Fraunhofer lines were removed from the Eris spectra by dividing 
each 600-second Eris spectrum by the spectrum of a solar analog star, SA 93—101. The 
airmass difference between each Eris spectrum and its corresponding solar analog spectrum 
was < 0.1. 

Since we saw no significant difference between the 24 600-second exposures of Eris, we 
summed the individual exposures to yield a single spectrum with a four hour exposure time. 
In Figure 1, we present the resulting reflectance spectrum. The spectrum was normalized 
to 1.0 between 8200 and 8400 A. The absorption bands in the spectrum are due to methane 
ice and the reddish slope in the "continuum" is consistent with absorption by a tholin-like 
material. 

To assess the uncertainty in our wavelength calibration, we measured the wavelengths 
of well- resolved sky emission lines in our spectrum a nd co mpared them to the VLT 



high— spectral resolution sky line atlas of 



Hanuschiki (120031 ). Unfortunately, due to the 



rather low spectral resolution of our spectrum, only the 6300.3 A and 9376.0 A sky lines 
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were well-resolved. Our measured centroids for these two lines were 0.4 and 0.8 A different 
from the centroids of Hanuschik. So, it appears our sky lines, and hence our wavelength 
measurements, are accurate to ~ 0.8 A. 

3. Laboratory Experiments 
3.1. Setup 

We performed experiments in the ice spectroscopy laboratory at Northern Arizona 
University. This facility includes a series of enclosed ice sample cells machined by the Lowell 
Observatory instrument shop out of aluminum alloy 7075 stock. Each cell was designed 
to be mounted upon the upward-pointing cold tip of a CTI Cryogenics model 1050C two 
stage closed cycle helium refrigerator driven by a CTI 9600 compressor. A 15 mm diameter 
cylindrical cavity in each cell holds the cryogenic ice sample for spectroscopic study. To 
permit transmission of the spectrometer beam through the samples, this cylindrical cavity is 
capped on each end with a sapphire window sealed in place with indium gaskets. Different 
cells have different spacing between their windows, enabling transmission spectroscopy 
through samples of various thicknesses. For the work presented here, we used a cell with a 
5 mm thickness. Sample materials were introduced and removed through a stainless steel 
flexible fill tube which entered the cell from above. 

Temperatures were measured using a pair of Lake Shore silicon diode thermometers, 
mounted on copper rods press-fitted into the cell. One of these is located just above the 
sample cavity and the other just below the sample cavity. A pair of 50 ohm nickel-chromium 
heater wires was wrapped around each cell and cemented into place with Emerson and 
Cuming's thermally conductive StyCast epoxy. As with the diodes, one heater was wrapped 
above and one was wrapped below the sample cavity. A Cryo-con model 44 temperature 
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controller was used to monitor and to control the temperature of the cell, via proportional, 
integral, and derivative (PID) closed-loop control of both heaters. By adjusting the relative 
amounts of power going into the upper and lower heaters, it was possible to induce a 
small vertical thermal gradient of up to approximately 1 °K/cm across the sample (2 °K 
between diode thermometers, with more power into the upper heater corresponded to a 
larger vertical temperature gradient). Temperature control stability with the Cryo-con 
44 was excellent, with unintended temperature fluctuations generally below 0.01 degree. 
There is a larger absolute uncertainty of less than 1 K on the sample temperature, owing 
to uncertainties in the calibration of the diode thermometers as well as the fact that they 
are mounted above and below the sample cavity, rather than directly sampling the ice. We 
were able to cool samples down to ~ 40°K. 

The cold head and cell were mounted within an evacuated and turbo-pumped stainless 
steel cylinder for insulation, and to eliminate contaminants which could condense on the 
exterior of the cell. The existence of the outer cylinder necessitates an additional pair of 
windows to pass the spectrometer beam. 

We used a Perkin Elmer Nexus 670 Fourier transform infrared (FTIR) spectrometer to 
generate an external beam, and glass lenses to focus the beam within the cell and again 
onto a silicon detector. All experiments covered a frequency (wavelength) range from 
~15000 cm" 1 (6667 A) to ~ 9000 cm" 1 (11111 A) at a spectral resolution of ~ 0.5 cm" 1 , 
although lower resolution was used in a few cases, where the signal to noise ratio was poor. 
At least 128 scans were taken in each experiment. In some experiments, 512 scans were 
taken to improve the signal to noise ratio. We present a schematic diagram of our lab set 
up in Figure 2. 



3.2. Sample Preparation 



Three separate gas bottles could be connected to a stainless-steel manifold through 
standard VCR vacuum fittings. In the experiments reported here, we used methane, 
nitrogen, and argon bottles. The purities of the gasses were 99.999% for methane and 99.9% 
for nitrogen and argon. A gas was released into the manifold, and then pumped out by a 
turbo pump and monitored by a 1000 Torr manometer. The purging process was repeated 
five to ten times. After the purges, the manifold was filled with a gas mixture dependent on 
the particular experiment. The cell was cooled to a temperature slightly above the melting 
point of the ice. In the case of mixtures, the temperature depended on the mixing ratio of 
the gases. The gases were then bled into the cell, at which time a liquid could be clearly 
seen to fill the cell. When the level of liquid was at least half the total volume of the cell, 
we stopped the gas flow. The manifold did not have a large ballast volume, and two fills 
of the manifold were usually required to achieve a satisfactory volume in the cell. After 
filling the cell with liquid, we slowly ramped the temperature down to below the melting 
point. During the ramp down, a gradient of 2 °K was typically kept on the cell. The target 
temperature was not always known, as the melting point for binary alloys varies a good 
deal. During the ramp down, the ice was periodically observed for freezing. If during the 
freezing process, the ice became optically thick, the process was stopped and modified until 
an optically clear ice could be made. Typically, we melted the optically thick ice and then 
reduced the size of the gradient across the sample and ramped the temperature down more 
slowly than the previous unsuccessful attempt to grow a clear ice. It should be noted that 
making an ice is trivial; however, making an optically clear, two-component ice is a difficult 
process. 

Spectra were taken during sample preparation as follows. A background spectrum was 
taken before we filled the cell with liquid. Another spectrum was taken after we filled the 
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cell with liquid. After growing a clear ice, we ramped the temperature down further, and 
spectra were taken at regular temperature intervals. 



3.3. Data Analysis 



Our spectra were analyzed as follows. Each spectrum was fitted with a LOWESS 



filter (j Cleveland 



19791 ) to take out some minor high frequency noise. We assumed that the 
absorption in our sample was governed by the Beer-Lambert Law, and so we computed 
Lambert absorption coefficient spectra, a(v). 

We measured the methane and nitrogen stoichiometry of a sample as follows. We 
integrated a band in a sample and the corresponding band in the pure methane sample 
to obtain integrated absorption coefficients. Next, we divided the integrated absorption 
coefficient of the sample by the integrated coefficient of pure methane. The result gave us 
the fraction of methane in our sample. We repeated the procedure for the other bands 
in the spectrum. From the scatter of the values for a sample, we found our solid-phase 
concentrations have an uncertainty of ~ 5% of the reported concentration. We found most 
of the uncertainty came from fitting backgrounds to the spectra. 

We computed the spectral shift for each band using a cross-correlation technique. We 
shifted a band in the pure methane spectrum relative to the same band in the sample 
spectrum over a range of shifts. We computed y 2 for each shift. We fitted a parabola to a 
plot of x 2 vs - shift, the minimum of the parabola gave us the shift for the band. Again, we 
found most of the uncertainty came from fitting background to the spectra. We found our 
shifts have an uncertainty of ~ 10% of the reported shift. 
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3.4. Methane-Nitrogen Results 



Next, we present our results for methane-nitrogen ice mixtures. In Figure 3, we plot 
the 11240 cm -1 (8897 A) band in the spectrum of a pure methane ice sample (black line). 
In addition, we plot the same band in the spectrum of a 6 % methane and a 94% nitrogen 
sample (red line). Both sa mples were at T = 60 °K. Our spectrum for pure methane is 



in excellent agreement with 



Grundy et al. 



(120021 ) . From cross correlation experiments, we 



found the 11240 cm -1 (8897 A) band in the spectrum of the 6% methane sample was blue 
shifted by 26 cm -1 (21 A) relative to the band of the pure methane sample. In Figure 4, 
we plot the absorption coefficient spectrum for the 13706 cm -1 (7296 A) band. We found 
the band in the spectrum of the 6% sample was blue shifted by 31 cm -1 (17 A) relative to 
the band of the pure methane sample. 

In Table 1, we present our blue shifts for samples with methane concentrations of 2, 
6, 24, and 99% relative to the bands of pure methane samples. In Figure 5, we plot the 
values in Table 1. Two trends are apparent in the table and figure. First, the larger the 
nitrogen concentration, (and so the smaller the methane concentration), the larger the blue 
shift. Second, t he higher frequency band always had a larger blue shift than the lower 



frequency band. 



Quirico fc Schmitt 



(119971 ) found the same trend of increasing blue shift 



Quirico fc Schmitt 



(1997) 



with increasing frequency in infrared bands. Finally, we note that 
found the 11240 cm -1 band had a blue shift of 33 cm -1 in spectra of samples with methane 
abundances of 0.1, 0.25, 0.8, and 2% at 44 °K. Their shift is in good agreement with our 
value of 29 cm -1 for a 2% methane sample (see Table 1). 

From Figure 5, it is clear that we did not study intermediate compositions. Methane 
rich gas mixtures did not easily go into the liquid phase. We had to dramatically increase 
the gas pressure and substantially lower the temperature, and then almost no liquid would 
condense, even at methane concentrations above 95 %. This effect was seen at all times, 
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and is mostly likely related to deviations from ideal mixtures. When we made a liquid, the 
ice grown was not optically clear. 

Besides studying the effect of concentration on shift, we wanted to perform an 
experiment to study the effect of temperature on shifts. Unfortunately, upon lowering the 
temperature the ice went through a phase transition (see Section 4.2), and it developed 
strong scattering characteristics. 

3.5. Methane- Argon Results 

Next, we present our results for methane- argon ice mixtures. In Figure 6, we plot the 
11240 cm" 1 (8897 A) band in the spectrum of a pure methane ice sample (same black line 
as in Figure 3). In addition, we plot the same band in the spectrum of a 9% methane and 
a 91% argon sample (red line). Both samples were at T = 60 °K. From cross correlation 
experiments, we found the 11240 cm -1 (8897 A) band in the spectrum of the 9% methane 
sample was blue shifted by 38 cm -1 (30 A) relative to the band of the pure methane sample. 

In Table 2, we present our blue shifts for bands in the spectra of samples with methane 
concentrations of 9, 15, 41, and 96% methane relative to the corresponding bands of pure 
methane samples. In Figure 7, we plot the values in Table 2. As in the methane-nitrogen 
samples, two trends are apparent in the figure and table. First, the larger the argon 
concentration, (and so the smaller the methane concentration), the larger the blue shift. 
Second, the higher frequency band always had a larger blue shift than the lower frequency 
band . 

Besides studying the effect of concentration on shifts, we performed an experiment 
to study the effect of temperature on shifts despite the ice developing strong scattering 
characteristics. Specifically, we cooled the 9% methane 91% argon ice from ~ 60 °K down 
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to about 40 °K. In Figure 6, we plot the 11240 cm" 1 (8897 A) band at T ~ 40 °K (blue 
line). It is clear that the single peak is split into two peaks. We warmed the sample up to 
T ~ 60 °K, and once again obtained a single peak (red line in Figure 6). We discuss the 
reason for the change in band profile with temperature in section 4.2. 



4. Models 

Next, we describe the models we used to derive methane and nitrogen abundances 
on Pluto and Eris. First, we describe a pure methane model. Then, we describe a binary- 
mixture model. 



4.1. Pure Methane 



In order to establish continuity with previous work, we compared our spectra of Pluto 
and Eris to spectra of pure methane ice at 30°K . We us ed Hapke theory to transform 



laboratory optical constants of pure methane ice at 30°K (iGrundy et al. 



20021 ) into spectra 



suitable for comparison to the Pluto and Eris spectra. In particular, Hapke t heory accoun ts 



for multiple scattering of light within a surface composed of particulate ice (IHapke 



19931 ). 



We used Hapke model parameters of h = 0.1, B G = 0.8, 9 = 30°, P(g) = a two component 
Henyey-Greenstein function with 80% in the forward scattering lobe and 20% in the back 
scattering lobe, and both lobes had asymmetry parameter a = 0.63. Our implementation 
of the Hapke model used two particle sizes and a tholin-like absorber. We explored 
parameters specific to our implementation, i.e. particle diameters (Di and D2), percent 
of large particles Di by volume (f), and tholin-like absorption (P x and P2), with Monte 
Carlo techniques (as described in more detail by Abernathy et al. 2009). For example, 
we performed cross-correlation experiments between 10,000 Hapke model spectra and the 
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portion of Pluto's spectrum containing the 8897 A methane band. We shifted each Hapke 
model spectrum in 2 A steps between —30 A and +0 A, calculated an adjusted x 2 goodness 
of fit, R, for each shift, plotted R vs. shift, and fitted a parabola to the plot to determine 
the shift corresponding to the minimum in R for each model. 

We found a dozen or so of 10,000 Hapke models gave equally good fits to Pluto's 8897 
A band. The relatively small differences in the parameters of these models had little effect 
on the blue shift. Our best fit model for the Pluto 8897 A band had grain sizes of Di=1.8 
cm (89% by volume) and D 2 =0.1 mm (11% by volume). It is important to recognize that 
our Hapke parameters do not represent a unique fit to Pluto's band; however, they are 
plausible values for transparent, pure metha ne ice grains and are comparable to Hapke 



parameters in previous fits of Pluto spectra (IGrundy fc Buie 



200lh . 



Our best fit model indicated Pluto's 8897 A band is blue shifted by 17 ± 2 A relative 
to pure methane ice. On 2007 June 2 1 UT , we were observing the anti-Charon facing 
hemisphere of Pluto. I Grundy &: Fink! (119961 ) reported a correlation between blue shift and 
longitude for Pluto's 8897 A band, and found a blue shift of ~ 12 A for the anti-Charon 
facing hemisphere. Grundy and Fink measured their shifts in a manner different than we 
described above and so differences between the two values may not be significant. 

In Figure 8, we present a comparison of Pluto's 8897 A band (black line) and our 
(shifted) best fit pure methane ice Hapke model (red line). It is clear that Pluto's 8897 A 
band is wider than the pure methane 8897 A band. We found \ 2 = 14. Unfortunately, 
noise in our Pluto spectrum prevented us from analyzing any other methane bands. 

We calculated the uncertainty in our blue shift as follows. The best fit (shifted) model 
was subtracted from Pluto's 8897 A band, giving us the noise in the astronomical spectrum. 
We note that the noise in the astronomical spectrum dominates noise in the Hapke model 
spectrum. Next, we calculated the standard deviation of the noise. Then, we applied 
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Gaussian noise with the same standard deviation to the model spectrum. Next, we applied 
our cross correlation technique to the noisy model and the original model spectrum. We 
generated 1000 different noisy model spectra and repeated the cross correlation experiment 
for each noisy model and the original model spectrum. A histogram of the resulting 1000 
shifts was fit with a Gaussian distribution, the standard deviation of the fit giving the 
uncertainty in our shift measurement, 2.0 A. 

We repeated the analysis described above on our Eris spectrum. In particular, we fit 
pure methane Hapke models to the 7296 A, 8691 A, and the 8897 A methane ice bands. In 
Table 3, we present our blue shift measurements for these bands. We point out that the 
Eris sh ifts in Table 3 are in g ood agreement with our previously reported shifts for these 



bands (Abernathy et al. 



20091 ) . In Figures 9 and 10, we present our Eris 8897 A and 7296 
A bands (black line) and our (shifted) best fit pure methane Hapke model (red line). We 
found x 2 = 10 for the 8897 A band and x 2 = 5 for the 7296 A band. 



4.2. Binary Mixture 

We began thinking about the mechanism responsible for the creation of the double 
peak in the 11240 cm -1 (8897 A) band upon cooling the methane-argon sample down 
(Figure 6), and it led us t o a n ew model. Examining the methane-argon phase diagram of 



Greer. Meyer, fc Barrett 



(119691 ) that we display in Figure 11, we found that cooling the 
sample down resulted in changing the structure of the ice from a soluble solution with 
methane and argon molecules mixed throughout the ice to a solution of limited solubility 
where we had highly methane rich regions and highly argon rich regions. The two peaks 
in the band came from the two compositionally distinct regions. Greer, Myer, & Barrett 
reported that at T ~ 40 °K, the solubility limit of methane in argon is ~ 13 % and the 
solubility limit of argon in methane is ~ 12 %. We believe this is the first time anyone has 
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spectroscopically observed the change in a methane-argon ice from a face-centered cubic 
(FCC) solid solution to regions dominated by FCC argon and FCC methane. 



The methane-nitrogen phase diagram of (jProkhvatilov fc Yantsevichlll983l ) that we 
display in Figure 12 is similar to the methane-argon phase diagram. At T ~ 60 °K our 
samples were solid solutions. At T ~ 40 °K, methane and nitrogen have limited solubility. 
The solubility limit of methane in nitrogen is ~ 2% and the solubility limit of nitrogen in 
methane is ~ 4%. Because of the experimental difficulties of growing an optically clear 
methane-nitrogen ice and cooling it, we were unable to spectroscopically see the structural 
change in a methane-nitrogen ice that we saw in the methane-argon ice. However, the 
methane-nitrogen phase diagram says it is there. 

At the temperatures of Pluto (T ~ 40°K) and Eris (T~30°K), we expect methane and 
nitrogen to have limited solubility in one another. So, we developed an ice model for Pluto 
and Eris consisting of two components — methane-rich crystals near the nitrogen solubility 
limit and ni t rogen -rich crystals near the methane solubility limit. We point out that 



DouteetaL 



(119991 ) made use of methane- rich regions and nitrogen-rich regions in modeling 
their near-infrared spectrum of Pluto. Although we could have added additional parameters 
to enable the two materials to have different grain sizes, we didn't feel the current data 
supported additional free parameters, and so we assumed the same grain sizes for both 
types of materials. In our model f is the overall methane abundance, r\ is the fraction of 
methane rich crystals, Sjv 2 is the solubility limit of nitrogen in methane, and Sch^ is the 
solubility limit of methane in nitrogen. The relationship between these quantities is given 
by 



/ = 7/(1 - 5jv 2 ) + (1 -y) S C h 4 ■ 
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Solving for 77, we get 



We combined absorption coefficients of the methane rich crystals (a u ), and absorption 
coefficients of the highly diluted methane crystals (a a ), 



= 77(1 - 5jv a )o! u (i/) + (1 - i])S C H 4 a s (u). (3) 

We approximated a u with pure methane coefficients, and a s with pure methane coefficients 
blue shifted by amou nts observed in highly di luted samples, e.g. 33 cm -1 for the 11240 



cm 1 (8897 A) band (IQuirico fc Schmitt 



19971 ). 



We fitted the Pluto and Eris bands much the same way we fitted them with pure 
methane models in the above section. The primary difference between the two models is 
that the binary mixture model had an additional free parameter, the methane abundance, 
f. In Figure 13, we plot our best fit binary mixture model (red line) to Pluto's 8897 A band 
(black line). Notice that the binary mixture model provides a better fit to Pluto's wide 8897 
A band than the pure methane model. We found x 2 = 10 f° r t ne binary mixture model 
and x 2 = 14 for the pure methane model. In addition, there was no apparent wavelength 
mismatch between the model and the data. The best fit model had a methane abundance 
of 3.2 ± 0.2 %, and hence a nitrogen abundance of 96.8%. 

In Figure 14, we plot our best fit binary mixture model to Eris' 8897 A band. Again, 
there was no apparent wavelength mismatch between the model and the data. The best fit 
model for Eris had a methane abundance of 10± 2 %, and hence a nitrogen abundance of 
90%. We found x 2 = 10 f° r t ne binary mixture model and x 2 = 10 for the pure methane 
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model. In Figure 15, we plot our best fit binary mixture model to the 7296 A band. We 
found x 2 = 3 for the binary mixture model and x 2 = 5 for the pure methane model. There 
was no wavelength mismatch between the model and the data. We found a methane 
abundance of 8 ± 2 %. In Table 4, we present the methane abundances corresponding to 
the best fit models for all three Eris bands. 

We analyzed our Eris bands with a methane-argon ice mixture. We found models with 
40 ± 5 % methane, and hence 60 % argon, gave the best fits to Eris' 7296 A, 8961 A, and 
8897 A bands. 



5. Discussion 



5.1. Eris' Stoichiometry 



Previous work suggested that the majority of methane ice on Eris is in a pure form 



rather than in a h i ghly diluted 



Dumas et al. 



2007 



Merlin et al 



brm w ith nitrogen ( IBrown. Trujillo. fc Rabinowitz 



2005 



20091 ) . Such thinking rests on two kinds of observations. 



First, the nitrogen 2.15//m band has yet to be detected in spectra of Eris, but it is seen in 
spectra of Pluto. Second, methane band shifts on Eris are smaller than on Pluto. 

Our analysis suggests nitrogen is almost an order of magnitude more abundant than 
methane on the surface of Eris. By fitting three methane bands in our Eris spectrum with 
the binary mixture model, we find a methane abundance of ~ 10 %, and hence a nitrogen 
abundance of ~ 90 %. Our model is able to fit all three bands. The lack of a 2.15/mi 
nitrogen band detection in spe ctra of Eris does not concern us too much. With an albedo 



of 70% ( iStansberry et al.l 120081 ) and a heliocentric distance of 97 AU, we expect a subsolar 
surface temperature of ~ 30°K for Eris. Such a temperature is below the T= 35.6 °K 
transition temperature between the a and (3 phase of nitrogen. So it is possible nitrogen is 
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in the a phase on Eris rather than the /3 phase as it is on Pl uto. The 2.15 \x absorption 



featu re in the a phase is much narrower than for the /3 phase (jGrundy. Schmitt. fc Quirico 



19931 ). If nitrogen is in the a phase, it will be more difficult to detect on Eris than on Pluto 
and it could co mplicate our analysi s. Spe cifically, the nitrogen in our experiments was in 



the (3 phase and 



Quirico Schmitt 



( 119971 ) found the v% + 1/4 methane band had shifts and 
profiles that strongly depended on the nitrogen phase whereas the V\ + v± methane band 
had shifts and profiles that were weakly dependent on the nitrogen phase. Our lab setup 
was not capable of going cold enough to grow nitrogen in the alpha phase and see whether 
any of our methane bands exhibit a strong nitrogen phase dependence. Another possibility 
is that some of the nitrogen ice may be hidden by the presence of scatterin g inclusions such 



Grundy et al. 



( 120101 ) suggested such 



as voids or wind blown dust particles of another ice. 
hidden nitrogen on Triton. Our analysis suggests that the volatile composition of Eris and 
Pluto may not be so different. 

Our spectra also are consistent with Eris having a surface composed of 60% Ar and 
40% CH4. It not unreasonable to expect A r on the surface of Eris, a lthough Ar is 30 times 



less abundant than N in the Solar System ( lAnders fc Grevesse 



1989J). If all of the N were 



in the form of N2, then Ar would be 15 times less abundant than N2. If some of the N 
were locked up in the interior of Eris in the form of NH3, the Ar abundance at the surface 
could be further boosted relative to N 2 . Furthermore, Ar is not destroyed by photolysis and 
radiolysis as is CH 4 . Eris could have started with a CH 4 dominated surface, but photolysis 
and radiolysis destroyed much of it, leaving Eris' surface dominated by Ar. Finally, the 
atomic mass of Ar is larger than the molecular masses of CH4 and N2, and so Eris would 
be able to retain Ar better than CH 4 and N 2 in its atmosphere and hence on its surface. 
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5.2. Eris Stratigraphy 



Previous work looked for a correlation between shift (a proxy for the nitrogen 
abundance) and albedo of maximum absorption (a proxy for depth into the surface) of 
methane bands . We and others reported such correlations, the signature of stratification. 



Licandro et al. 



(120061 ) found essentially no shift for the weaker 7296 A band, and a shift 
of 15 A for the stronger 8897 A band. Such a trend suggested a decrease in the nitrogen 
abundance with depth into the surface. On the other hand, we reported an increase in 



nitrogen abundance with depth ( lAbernathy et al 



20091 ). In the work we report here, we fit 



models with the same methane and nitrogen stoichiometry for three bands. Furthermore, if 
there were a change in the stoichiometry with depth, the weaker bands which probe deeper 
into the surface and therefore probe a range of mixing ratios, should be broader than the 
stronger bands which don't probe as deep into the surface. We find no evidence that the 
weaker 7296 A and 8961 A bands are systematically broader than the 8897 A band. In 
short, we see no evidence of stratification in the ice. 



5.3. Pluto's Stoichiometry 



There are a range of reported values for the methane and nitrogen stoichiometry of 



Pluto in the literature. 



Owen et al. 



(119931 ) reported a methane abundance of 1.5% and an 



uncerta inty of a factor of t wo. So the Owen et al. methane abundanc e could be as 



a 3 %. iDoute et al. 



( 119991 ) found a methane abundance of 0.5%, and 



Olkin et al 



arge as 



(120070 



found a methane abundance of 0.36%. 



Our methane abundance of 3.2 ± 0.2 % overlaps with the values reported by 



Owen et al. 



( 119931 ). However, our methane abundance is inconsistent with abundances smaller than 



2%. At such low abundances, the ice would be dominated by highly diluted methane 
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crystals. So, we should see an 8897 A band that is as narrow as the pure methane 8897 A 
band, but shifted 26 A (33 cm -1 ) relative to the pure methane band. Clearly, the Pluto 
band in Figure 8 is much wider than the pure methane band and not shifted enough for 
highly diluted methane. We point out that the lOlkin et al.l ( 120071 ) methane abundance of 



0.36% is based on 3 jum data, and so does not probe as deep into the surface as our optical 
data. Perhaps a change in composition with depth is responsible for the different values. 



5.4. Future Work 

We hope to obtain near- and mid-infrared spectra of Eris at higher spectral resolutions 
than reported in the literature (R > 15,000). Such observations will test the stoichiometry 
and the lack of stratification in Eris' ice described here. In addition, we hope to obtain 
additional spectra of Pluto and map its methane and nitrogen abundances as a function 
of longitude. Such measurements may prove useful in the interpretation of New Horizons 
observations of Pluto. 
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Wavelength (Angstroms) 

Fig. I. — Spectrum of Pluto (top) taken on 2007 June 21 UT with the Steward Observatory 
2.3 m telescope, normalized to 1 between 8200 A and 8400 A and shifted upward by 0.2 
to facilitate a comparison with the spectrum of Eris. Spectrum of Eris (bottom) taken on 
2008 October 3 UT with the 6.5 m MMT, and normalized to 1 between 8200 A and 8400 
A. The absorption bands in both spectra are due to methane ice and the reddish slope in 
the "continuum" is consistent with absorption by a tholin-like material. Eris has deeper 
methane bands than Pluto. 
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Fig. 2. — A schematic diagram of the set up in our ice laboratory. 1 gas bottles; 2 turbo 
pump for vacuum that surrounds cell; 3 closed cycle helium refrigerator and compressor; 
4 temperature controller and monitor for cell; 5 Fourier transform spectrometer; 6 silicon 
detector; 7 to data acquisition computer; 8 cell. 
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Fig. 3. — The 11240 cm -1 (8897 A) band in the spectrum of a pure methane ice sample 
(black line). The same band in the spectrum of a sample of 6% methane and 94% nitrogen 
(red line). The band corresponding to the highly diluted methane sample (red line) has the 
same shape as the band corresponding to the pure sample (black line), but it is blue shifted 
by 26 cm" 1 (21 A). 
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Fig. 4. — The 13706 cm^ 1 (7296 A) band in the spectrum of a pure methane ice sample 
(black line). The same band in the spectrum of a sample of 6% methane and 94% nitrogen 
(red line). The band corresponding to the highly diluted methane sample (red line) has the 
same shape as the band corresponding to the pure sample (black line), but it is blue shifted 
by 31 cm -1 (17 A). The low signal to noise ratio of the diluted band is due to the weak 
absorption efficiency of the band and the low methane abundance of the sample. 
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Fig. 5.— Blue shift of the 11240 cm" 1 (8897 A) band (open circles) and the 13706 cm" 1 
(7296 A) band (filled squares) vs. methane abundance. The larger the nitrogen abundance 
(smaller the methane abundance), the larger the blue shift of the band. For the same 
methane concentration, the 13706 cur 1 band has a larger shift than the 11240 cur 1 band. 
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Fig. 6. — The 11240 cm^ 1 (8897 A) band in the spectrum of a pure methane ice sample at 
T ~ 60 °K (black line). The same band in the spectrum of a sample of 9% methane and 
91% argon at T ~ 60 °K (red line). The highly diluted methane band, but at T ~ 40 °K 
(blue line). At the higher temperature, the methane and argon molecules are soluble and 
methane and argon molecules are mixed throughout the ice. At the lower temperature, the 
methane and argon are not soluble and there are regions of nearly pure methane and regions 
of highly dilute methane giving rise to the two peaks in the band. 
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Fig. 7.— Blue shift of the 11240 cm" 1 (8897 A) band (open circles) and the 13706 cm" 1 
(7296 A) band (filled squares) vs. methane abundance. The larger the argon abundance 
(smaller the methane abundance), the larger the blue shift of the band. For the same 
methane concentration, the 13706 cur 1 band has a larger shift than the 11240 cur 1 band. 
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Fig. 8. — Pluto's 8897 A methane band taken on 2007 June 21 UT with the Steward 
Observatory 2.3 m telescope (black line), and the best fit pure methane Hapke model (red 
line). Blue shifting the pure methane band by 17 ± 2 A gave the best fit to Pluto's band. 
Notice that the Pluto band has a wider profile than the pure methane band, x 2 = 14- 
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Fig. 9.— Eris' 8897 A methane band taken on 2008 October 3 UT with the MMT 6.5 m 
telescope (black line), and the best fit pure methane Hapke model (red line). Blue shifting 
the pure methane band by 4.0 ± 0.6 A gave the best fit to Eris' band, \ 2 — 10 




Fig. 10.— Eris' 7296 A methane band taken on 2008 October 3 UT with the MMT 6.5 m 
telescope (black line), and the best fit pure methane Hapke model (red line). Blue shifting 
the pure methane band by 4.8 ± 1.0 A gave the best fit to Eris' band, x 2 = 5 
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Fig. 11. — Phase diagram for methane- argon ice mixtures ( IGreer. Meyer, fc Barrett 
In our experiments, we cooled an ice sample with ~ 91% argon and ~ 9% methane from 
60 °K to 40 °K. We believe the sample went from a FCC solid solution with methane and 
argon mixed throughout the sample to a FCC argon + FCC methane structure composed 
of highly methane rich regions and highly argon rich regions. 
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Fig. 12. — Phase diagram for methane-nitrogen ice mixtures ( jProkhvatilov fc Yantsevich 
19831 ). At the temperatures of our experiments (T ~ 60 °K), the samples were a solid 
solution of methane and nitrogen molecules (CH 4 I and /3 N 2 ). At the temperature of Pluto 
(T ~ 40 °K), the ice probably is composed of methane rich regions and nitrogen rich regions 
(CH 4 I + (3 N 2 ). At the temperature of Eris (T ~ 30 °K), the ice also probably is composed 



of methane rich regions and nitrogen rich regions (CH 4 I + a N 2 ). 
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Fig. 13. — Pluto's 8897 A methane band taken on 2007 June 21 UT with the Steward 
Observatory 2.3 m telescope (black line), and the best fit binary mixture Hapke model (red 
line). The best fit model had a methane abundance of ~ 3% and a nitrogen abundance of 
~ 97%. x 2 =10. 
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Fig. 14.— Eris' 8897 A methane band taken on 2008 October 3 UT with the MMT 6.5 m 
telescope (black line), and the best fit binary mixture Hapke model (red line). The best fit 
model had a methane abundance of 10% and a nitrogen abundance of 90%. x 2 = 10. 
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Fig. 15.— Eris' 7296 A methane band taken on 2008 October 3 UT with the MMT 6.5 m 
telescope (black line), and the best fit binary mixture Hapke model (red line). The best fit 
model had a methane abundance of 8% and a nitrogen abundance of 92%. \ 2 — 3- 
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Table 1. Band Shifts For Methane-Nitrogen Samples" 



Methane 6 


13706 c 


11240 c 


2 


34 


29 


6 


31 


26 


24 


25 


19 


99 


1 






a Samples at T = 60°K 
b Percent Methane 
c Shifts in units of cm" 
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Table 2. Band Shifts For Methane- Argon Samples" 



Methane 6 


13706 c 


11240 c 


9 


56 


38 


15 


46 


33 


41 


27 


20 


96 


1 


1 



a Samples at T = 60°K 
b Percent Methane 
c Shifts in units of cm" 
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Table 3. Band Shifts Relative to Pure Methane 



Object 7296 a 8691 a 8897' 



Pluto 17 ± 2 

Eris 4.8 ± 1.0 4.8 ± 1.1 4.0 ± 0.6 



a Shifts in units of A 
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Table 4. Methane Abundances" 



Object 7296 8691 8897 



Pluto 3.2 ± 0.2 

Eris 8 ± 2 9±2 10 ±2 



a Percent Methane 



